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Whereas  the  USSR  models  jet  or  plate  collapse,  under  shock  loading  conditions, 
via  visco-plastic  or  rate-dependent  models. 

The  USA  establishes  a jet-no-jet  criterion  and  a jet  cohesiveness  criterion 
based  on  Mach  number  considerations  or  impinging  angle  considerations.  Usually, 
the  Mach  number  is  based  on  an  ambient  bulk  speed  of  sound  value.  The  USSR, with 
their  incompressible,  viscid  flow  model,  uses  a critical  Reynolds  number 
criterion  to  define  the  jet-no-jet  condition.  7o  ui-e 

To  utilize  ^the  USSR  visco-plastic  model,  the  constant  (but  unknown)  viscosity 
coefficient  must  be  determined.  This  parameter  has  been  obtained  from  experi- 
mental data  for  various  materials  under  shock  loading  conditions.  The  purpose 
of  this  report  is  to  present  the  one-dimensional,  USSR,  viscous  model  for  jet 
or  plate  collapse  and  the  USSR  criterion  for  jet  formation  and  jet  cohesiveness. 
Also,  viscosity  values  for  various  materials  under  shock  loading  conditions  are 
given  from  the  USSR  open  literature  and  from  recent  USA  publications. 
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I.  INTRODUCTION 

The  USSR  has  been  engaged  in  the  study  of  visco-plastic  modeling  of  t 
shock  loaded  materials  since  1940  with  the  pioneering  work  of  Il'yushin1  ,2 
and  Popov3.  The  USSR  developed  models  to  describe  the  jet  formation 
process  during  the  collisions  of  metal  plates1*-9.  These  models  are 
analogous  to  those  employed  by  the  BRL10-12,  but  the  USSR  has  extended 
certain  models  to  include  the  impact  of  asymmetric  plates5  and  viscous - 
flow  effects1**  5.  The  basic  BRL  models10-12  employ  a one -dimensional, 
incompressible,  inviscid  flow  model  for  the  jet  or  plate  collapse  process, 
whereas  the  USSR  models4'5  treat  a one-dimensional,  incompressible,  vis- 
cous flow  model  for  the  shock  loading  of  metal  plates.  As  a result,  the 
criterion  used  for  predicting  a jet-no-jet  or  jet  cohesiveness  condition 
is  based  on  a critical  Mach  number  or  impingement  angle  criterion13'14. 
Often  the  critical  Mach  number  is  calculated  from  the  ambient  bulk  speed 
of  sound  value14.  The  USSR  models,  however,  use  a critical  Reynolds 
number  criterion  to  establish  the  jet-no-jet  condition4'5.  The  USSR 
viscous  jet  collapse  model  is  presented  in  Section  II. 

In  the  axisymmetric  hydrocode  models  used  at  BRL,  compressible, 
inviscid  flow  is  assumed,  but  the  constitutive  relationships  are  based  on 
elastic-perfectly  plastic  or  on  work  hardening  models.  The  USSR  utilizes 
rate-dependen'c,  visco-plastic,  stress-strain  relationships3-5,15-17. 

These  relationships  require  a knowledge  of  the  dynamic  viscosity  coeffi- 
cient, and  many  USSR  investigators1-4' 18-25  have  deduced  the  viscosity 
from  experimental  measurements  under  shock  loading  conditions.  Some 
Western  World  viscosity  values  have  also  been  obtained26-30  from  shock 
loading  experiments. 

The  dynamic  viscosity  is  dependent  on  many  parameters,  primarily 
strain  rate,  pressure,  and  temperature  and  can  range  from  10-3  Pa-s  for 
metals  in  the  liquid  state26  to  107  Pa-s  for  very  low  strain  rates4. 
However,  for  a given  pressure,  temperature,  and  strain  rate,  all  materials 
may  have  nearly  the  same  viscosity18  and  behave  as  Newtonian  fluids  (for 
a fixed  strain  rate)1-4'18-23.  A summary  of  the  experimental  and  theoreti- 
cal values  of  the  viscosity  coefficients  are  given  in  Section  III.  The 
references  cited  relate  to  viscosity  data,  visco-plastic  modeling,  shock- 
wave  propagation  and  dislocation  dynamics.  The  cited  references  repre- 
sent a partial  bibliography  which  is  far  from  complete. 


II.  VISCOUS  JET  FORMATION  MODEL 

Godunov,  Deribas,  and  Mali4  modified  the  jet  formation  equations  of 
References  10-12  to  include  the  viscosity  of  metals  under  symmetrical 
dynamic  loads.  The  resulting  equations  are  analogous  to  the  inviscid 
models10-12  and  include  extra  terms  resulting  from  viscous  flow  considera- 
tions. All  other  assumptions  are  the  same  as  employed  by  Reference  10. 


*Refevenoes  are  listed  on  -page  35 


The  equations  from  Reference  4 are  repeated  here  since  this  author  feels 
that  some  of  the  results  given  in  Reference  4 are  in  error.  Also,  an 
outline  of  the  derivation  and  the  analogy  to  References  10-12  are 
presented. 

Birkhoff10  gives 


cos  ~ (6 -a)  cos  j (B-a)  , 

V = V { ■ =— 3 + -i— r + sin  tc  (6-a) 

sm  B tan  B 2 K J 


(1) 


and 


n cos  ~ CB-ot) 

= V 


cos  a o sin  (B-a) 


(2) 


where  the  notation  of  Reference  4 will  be  converted  to  the  notation  used 
in  Reference  10,  and 


and 


V is  the  forward  jet  velocity, 

Vq  is  the  velocity  of  the  liner, 

D is  the  speed  of  the  plane  detonation  wave 
traveling  parallel  to  the  jet  axis, 

2a  is  the  initial  conical  apex  angle, 

2 B is  the  collison  angle. 


D1  is  defined  to  be  the  speed  of  the  detonation  wave  traveling  along 
the  liner  or 


d'=  ° 


cos  a 

The  combination  of  Equations  1,  2,  and  3 yields 


V = D sin  (B-a) 


1 + cos  B * 1 s 

— + tan  2 (S'a) 


(3) 


(4) 


1* 

This  result  disagrees  with  Equation  1 of  Godunov 4 apparently  due  to  an 
error  by  Godunov \ 
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The  mass  per  unit  length  of  the  jet  is  given  by 

».  - (s) 

3 ^ 

where  m is  the  mass  per  unit  length  of  the  collapsing  liner.  The  jet 
mass  equation  is  identical  to  that  given  by  Birkhoff10.  Also,  Godunov4 
gives,  for  the  jet  velocity, 


V. 

3 


v (1+cos  a) 
o sin  a 


(6) 


when  the  detonation  wave  moves  normal  to  the  surface  of  the  liner  or 
6 = a . This  equation  was  given  by  Birkhoff10. 

For  the  case  of  a liner  being  collapsed  by  a detonation  wave  moving 
along  its  surface, 


V = V + V*  = D s*n 
j c sin  8 


1 + 


cos 


cos 


( 

( 


~ 2 } 


(7) 


where  V is  the  velocity  of  the  stagnation  point,  and  V1  is  the  liner 
flow  veJocity  or 


V 

c 


sin  (B-ot) 
sin  8 


(8) 


and 


v'  - d sin  8 - sin  a (9) 

sin  6 

Equations  8 and  9 follow  from  De  Foumeaux's12  Equations  52  and  53, 
respectively,  for  the  stationary  or  steady-state  case  where  8 = a + <j> 
and  <j>  is  the  plate  bending  angle. 

Godunov4  proposed  a method  to  determine  the  influence  of  viscosity 
on  the  jet  formation  process.  In  particular,  a jet  formation  criterion 
is  derived  which  is  based  not  on  taking  account  of  the  compressibility, 
as  in  References  31  or  13,  or  on  taking  account  a critical  Mach  number, 
as  in  Reference  14,  but  on  taking  into  account  the  viscous  properties 
of  metals.  An  approximate  method  is  used  to  estimate  the  effect  of 
viscosity  in  the  plane  problem  of  jet  collisions.  The  fluid  is  incom- 
pressible, the  motion  irrotational  and  steady-state,  and  the  coefficient 
of  viscosity  is  constant.  Then  the  solutions  of  the  Euler  equations 
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automatically  satisfy  the  Navier-Stokes  equations  of  motion  and  the  whole 
difference  between  the  problems  of  ideal  and  viscous  jet  collisions  is 
the  conditions  on  the  free  surface  of  the  jet.  For  the  flow  fields  to 
agree  in  these  problems,  some  forces  must  be  applied  to  the  free  ,urface 
in  the  viscous  flow  case.  The  taking  into  account  of  the  influence  of 
these  forces  on  the  flow  from  collisions  between  jets  of  an  ideal  fluid 
will  be  the  estimate,  in  a first  approximation,  of  the  influence  of 
viscosity  in  the  jet  formation  problem4. 

The  components  of  the  viscous  stress  tensor  in  an  incompressible 
fluid  are  defined  by 


where  are  the  velocity  components,  x.  are  the  Cartesian  coordinates, 
and  u is  the  viscosity  coefficient.  The  modulus  of  the  viscous  force  is 
taken  to  be 

ui 

0 + i a | = 2y  — , (11) 

1 x y1  R 

where  u^  is  the  modulus  of  the  complex  velocity,  ui-iug,  and  R is  the 
radius  of  curvature  of  the  free  surface  of  the  fluid.  Then,  Reference 
4 postulates  that,  if  the  decelerating  effect  of  the  surface  viscous 
forces  on  the  reverse  jet  (or  slug)  will  equal  the  computed  value  of  the 
reverse  jet  force,  then  it  is  natural  to  expect  that  the  reverse  jet 
will  not  be  formed.  The  horizontal  component  of  the  viscous  force  acting 
on  the  free  surface  of  the  reverse  jet  over  the  angular  variation  from 
8 to  tt/2  is 


M1  = 2u  Uj  (1-sin  8) , (12) 

and,  since  the  same  force  per  unit  length  acts  on  the  opposite  symmetric 
part  of  the  free  surface  of  the  reverse  jet,  the  total  deceleration 
force  per  unit  length  is 

M = 2 M1  = 4p  Uj  (1-sin  6) . (13) 

Now,  consider  the  symmetric  collision  of  two  plane  fluid  jets.  For 
an  ideal,  incompressible  fluid  the  reverse  jet  force  per  un:t  length  is 
given  by  Reference  32  as 

Mj  = 2p  (5j  Uj2  sin2  8,  (14) 

where  p is  the  density,  5.  is  the  thickness  of  the  colliding  jets,  Uj 
is  the  reverse  jet  velocity  in  a coordinate  system  coupled  to  the 
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stagnation  point,  and  28  is  the  collision  angle.  Godunov4  assumes  that 
viscous  fluid  jets  collide  and  form  a reverse  jet  and  the  viscous  fluid 
flow  domain  coincides  with  the  analogous  domain  for  an  ideal  fluid.  As 
mentioned  previously,  the  solutions  of  the  Euler  equations  satisfy  the 
Navier-Stokes  equations  of  motion;  the  differences  in  these  flows  are 
attributed  to  the  conditions  on  the  free  surface. 


The  formation  of  a reverse  jet  is  possible  only  if  the  horizontal 
component  of  the  viscous  force,  acting  on  the  free  surfaces  of  the 
reverse  jet,  is  less  than  the  force  resulting  from  the  symmetric  collision 
of  two  plane  fluid  jets.  From  Equations  13  and  14,  this  inequality  is 
established: 


2 2 

2p  5^  u^  sin  8 > 4y  u^  (1-sin  8) . 
By  letting  the  Reynolds  number  be  defined  as 


(15) 


Re  = 


2 

6 u^  sin  8 
v (1-sin  8) 


(16) 


where  v = y/p  is  the  kinematic  viscosity,  the  criterion  for  jet  formation 
is 

Re  > 2.  (17) 

Next,  the  jet  velocity  for  a viscous  flow  is  computed  for  the  case 
where  a reverse  jet  exists.  A simple  force  balance  yields 

2p  u ^ sin2  8 = 2p  6^  u^'  sin2  8 - 4y  Uj  (1-sin  8),  (18) 

where  ^ now  represents  the  reverse  jet  velocity  in  a coordinate  system 
coupled  to  the  stagnation  point  and  uj  is  the  liner  flow  velocity.  From 
Equation  18  and  Equation  9,  the  reverse  jet  velocity,  U2>  for  the  impact 
of  two  flat  plates  is  given  by 


u ■ D 


(sin  8-sin  a\  [7~ 
sin  8 /Y  " 


2v  (1-sin  8) 

D sin  8 (sin  8-sin  a) 


(19) 


and  the  jet  velocity  in  a fixed  (laboratory)  coordinate  system  is 
obtained  by  adding  V from  Equation  8 to  Equation  19,  or  for  D 
perpendicular  to  thecliner, 


13 


sin  p 


/•P-CU 
COS  (p-) 


Oj  D sm  B (sm  B-sin  a) 


For  plates  in  parallel,  or  when  a = 0,  Equation  20  becomes 


u.  = D 
3 


1 + - 1 - 


2v  (1-sin  g) 
2 

6^  D sin  g 


Equation  20  (or  Equation  21,  if  a = 0)  gives  the  shaped-charge  jet  velocity 
for  viscous  flow  effects. 

( When  the  detonation  wave  velocity  is  directed  along  the  liner,  i.e., 

D = D/cos  a from  Equation  3,  Equation  20  becomes  , for  D*  directed  alone 
the  jet  s 


u = D'  sin  (B-a)  , 003  ( 2 } 

j cos  a sin  B ,g-cu 

cos  (-p) 


L 2v  (1-sin  B)  cos  a I 

<$,  D'  sin  B (sin  B-sin  a)  j 5 


which  is  identical  to  Equation  21  when  a = 0. 

For  the  detonation  wave  impacting  normal  to  the  liner  surface,  the 
jet  velocity  from  Equation  22  for  B = ct  and,  by  using  Equation  2, 
becomes 


u.  - A_  [l  + cos  a Jl  - 
J sm  a W 6.  V sin  2 a 


(23) 


Now,  Equations  23,  22,  and  20  relax  to  the  inviscid  jet  relationships 
given  by  Equations  6,  4,  and  7,  respectively,  for  the  detonation  wave 
directed  along  the  liner  when  v = 0. 

The  values  of  the- metal  viscosity  coefficients  in  shaped-charge  jet 
flows  must  be  used  in  conjunction  with  the  strain  rates  achievable  which 
can  be  estimated  as 


. u.  (ir-3) 

• Au  1 v J 

As  'R-(ir-B) 


2 

2 6^  sin  3 


(24) 


where  a change  in  velocity  results  from  the  rotation  of  the  velocity 
vector,  Uj,  through  the  angle,  ir-g,  over  the  path.  As  = R (ir-B).  The 
radius  of 1 curvature,  R,  is  given  as 


R = 


2 6 


1 • 2 . 
- sm  6 


(25) 


from  Reference  9 where  the  3 given  in  Reference  9 corresponds  to  23  in 
the  notation  of  Reference  4. 


Final: 7 . experimental  jet  data  can  be  used  to  estimate  the  steady 
state  visc.-sity  coefficient,  where,  from  Equation  18, 


y 


„ . 2 
p u^  sin  3 

2 (1-sin  3) 


(26) 


Also,  from  Equation  15,  the  critical  velocity  of  the  liner  in  the 
direction  of  the  stagnation  point,  at  which  a failure  to  jet  will  occur, 
is  given  by 


cr  . 2 v . (27) 

1 61  sin2  3 

Finally  Equations  17,  20,  26,  and  27  define  the  USSR  criterion  for  jet 
formation  and  give  the  shaped-charge  jet  velocity  for  viscous  effects. 
Again,  Godunov^  was  the  major  source  of  this  information.  This  model  has 
been  extended  to  asymmetric  collisions  in  Reference  5.  Reference  6 
derives  an  expression  for  the  stagnation  point  velocity,  in  steady  state, 
as  a function  of  the  charge-to-mass  ratio  for  inviscid  flow. 
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III.  VISCOSITY  DATA 

The  earliest  known  USSR  viscosity  measurements  were  from  Il'yushin1*2 
and  Popov3  as  reported  in  Reference  18.  The  material  behavior  is  modeled 
as  a rate-dependent  process  and  the  visco-plastic  relationship  between 
the  stress,  a,  and  the  strain  rate,  e,  is  given  by 

a = a + p e, 

y 

where  Oy  is  the  yield  stress  and  p is  the  dynamic  viscosity  coefficient. 
Thus,  knowledge  of  the  viscosity  coefficient,  which  is  assumed  to  be 
constant,  is  required  for  the  application  of  the  visco-plastic,  stress- 
strain  rate  relationship.  As  will  been  seen,  the  viscosity  coefficient  is 
in  reality,  not  constant  but  dependent  primarily  on  the  pressure, 
temperature,  and  strain  rate.  Nonetheless,  Il'yushin1*2  and  Popov3 
performed  experiments  on  cylindrical  specimens  tested  in  a pneumatic 
impact  tester  and  concluded  that  the  viscosity  coefficients  of  various 
steels  lie  in  the  interval  of  3 x 104  to  14  x 104  Pa-s  and  the  viscosity 
coefficients  of  various  aluminums  lie  in  the  interval  of  3 x 104  to 
4 x 104  Pa-s. 

Another  method  of  determining  the  viscosity  of  continuous  media, 
based  on  the  experimental  investigation  of  the  development  of  small 
perturbations  in  a shock  front,  was  proposed  ir.  Reference  23. 

Sakharov23  attempted  to  investigate  the  stability  of  a plane  shock 
wave  in  a substance  in  a condensed  state  and  to  deduce  the  mechanical 
properties  of  the  substance  under  conditions  of  high  pressure  and  tem- 
perature behind  the  shock  front.  The  experimental  arrangement  is 
depicted  in  Figure  1 from  Reference  23.  The  experiment  was  designed  to 
satisfy  the  boundary  and  initial  conditions  associated  with  theoretical 
calculations  performed  for  a linear  approximation  (ka  <<  1 where  a is 
the  perturbation  amplitude,  k = 2tt/A  and  X is  the  wavelength).  For 
t = t0,  the  surface  of  the  shock  front  is  assumed  to  have  a sinusoidal 
profile  and  the  flow  behind  and  ahead  of  the  shock  front  is  constant. 

From  Figure  1,  a plane  shock  wave,  originating  from  the  detonation 
of  the  explosive  charge,  passes  thi'ough  the  sinusoidal  grooves  of  the 
disc  (2)  into  the  wedge  (3) . Perturbations  appear  on  the  wedge  of  the 
same  wavelength  as  the  grooves  in  the  disc.  The  wedge  and  the  disc  are 
made  of  the  material  under  investigation.  On  further  propagation  of  the 
shock  wave  through  the  wedge  which  is  carrying  the  perturbation,  the 
shock  wave  enters  the  gap  between  the  tapered  surface  of  the  wedge  and 
the  plastic  plate,  where  the  perturbations  are  recorded  by  a fluoro- 
metric  device  SFR-2M.  By  using  slits  which  are  perpendicular  to  the 
direction  of  the  wedge  taper,  a single  perturbation  can  be  recorded  in 
one  experiment  at  several  successive  times23. 
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By  choosing  sufficiently  large  values  for  the  length  and  diameter 
of  the  explosive  charge,  the  diameter  of  the  disc,  etc.,  it  is  possible 
to  approximately  insure  a uniform  flow  behind  the  shock  front23. 

The  method  of  producing  the  perturbations  destroys  the  uniformity  of 
flow  behind  the  shock  front  at  distances  of  order  aQ,  the  initial  ampli- 
tude, which  leads  to  errors  of  the  second  order  with  respect  to  kaQ.  In 
each  case,  the  effect  of  deviations  from  the  calculated  conditions  in  the 
experimental  arrangements  was  studied  experimentally,  and  the  experimental 
results  were  subjected  to  corrections23. 

Sakharov23  studied  the  development  of  perturbations  by  shock  waves  in 
aluminum  alloy  AL-9  (90%  aluminum).  The  pressure  at  the  wave  front  was 
3.14  x 1010  Pa,  the  temperature  was  603°K,and  the  density  was  3.4  Mg/m3. 
These  parameters  were  obtained  from  an  assumed  equation  of  state33. 

The  experimental  curves  of  the  development  of  perturbations  were  obtained 
for  wavelengths  of  0.1,  0.2,  and  0.33mm  and  initial  perturbation 
amplitudes  of  kaD  = 0.29  to  1.74.  In  the  experiments  using  different 
wavelengths,  all  linear  parameters  of  the  experimental  assembly  were 
scaled  with  the  wavelength  to  eliminate  the  effect  of  boundary  conditions 
on  the  required  dependence  of  X.  The  viscosity  coefficient  for  the  AL-9 
aluminum  alloy  was  calculated  to  be  2 x 103  Pa-s  behind  the  shock  front23. 
Increasing  the  shock  wave  pressure  up  to  1.01  x 1011Pa  changes  the  coef- 
ficient of  viscosity  weakly,  but  it  does  not  exceed  104  Pa-s23.  The 
perturbation  development  curves  for  different  values  of  the  wavelength 
agreed  for  the  case  of  aluminum  powder  with  an  initial  density  of  0.68 
Mg/m3.  It  is  possible  that  this  indicates  melting  of  the  aluminum 
behind  the  shock  front,  since  the  use  of  a powdered  material  leads  to 
considerable  heating  behind  the  shock  front23. 

Later  Mineev21  continued  the  experiments  in  order  to  obtain  the 
viscosity  of  aluminum,  lead,  and  sodium  chloride  at  pressures  of  10^  to 
2 x 10i:iPa.  It  was  shown  that,  when  complete  geometric  similarity  is 
maintained  in  the  experiments  (with  respect  to  the  wavelength  of  the 
disturbances),  the  influence  of  viscosity  reduces  to  a phase  shift  be- 
tween the  curves  of  disturbance  amplitude  versus  time,  thereby  showing 
the  development  of  disturbances  having  various  wavelengths.  However, 
under  certain  thermodynamic  conditions,  the  development  of  a disturbance 
is  independent  of  the  wavelength.  This  is  associated  with  the  melting 
of  the  test  material  behind  the  shock  front. 


In  the  actual  experiments,  the  development  of  disturbances  of  two 
wavelengths,  namely,  Xj  = ,01mm  and  X2  15  .02mm  with  relative  initial 
amplitudes  of  ka^  =0.63  and  1.89,  respectively,  was  considered  in  two 
geometrically  similar  arrangements  with  respect  to  X. 
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The  dynamic  viscosity,  y,  was  obtained  from  the  formula"*- 
P = PD  Ax/kaQ  (l/Aj-l/A^, 

where  p is  the  density  of  the  material  behind  the  unperturbed  shock 
front,  D is  the  wave  speed  of  the  unperturbed  shock  front,  and  x = 

S(t)/A,  where  S is  the  path  traveled  by  the  shock  wave.  In  fact,  Ax  is 
the  phase  shift  of  the  curves  showing  the  development  of  disturbances 
with  Xi  and  A2,  plotted  with  cooidinates  y = a(t)/a0  and  x.  Here, 
a(t)  and  aQ  are  the  running  and  initial  disturbance  amplitudes,  respec- 
tively21. 

Table  I,  from  Reference  21,  lists  the  values  of  the  dynamic  viscosity 
for  lead  and  aluminum  under  various  pressures,  compressions,  and  tempera- 
tures. The  viscosity  values  are  given  with  a maximum  error  band  due  to 
the  inaccuracy  in  the  determination  of  Ax.  Again,  equations  of  state 
were  used  to  determine  the  thermodynamic  properties33"35. 

Also,  disturbances  in  the  shock  front  in  sodium  chloride  were  g 
investigated  for  two  cases  corresponding  to  a pressure  of  240  x 10°  Pa, 
a compression  of  1.44,  and  a temperature  of  1550°K,  and  a pressure  of 
205  x 108  Pa,  a compression  of  1.31,  and  a temperature  of  2620°K.  Again, 
empirical  equations  of  state  were  required36.  In  the  first  case,  the 
viscosity  was  deduced  to  be  2 ± 1 x 104  Pa-s  and  less  than  103  Pa-s  in 
the  second  case. 

From  the  data  presented  in  Table  I,  Mineev21  states  that  the  viscosity 
values  of  aluminum  and  lead  behind  the  shock  front  at  pressures  of  310 
and  350  x 108  Pa  are  nearly  equal  at  2-4  x 103  Pa-s.  Under  these  con- 
ditions, the  aluminum  and  lead  are  deformed  by  174%  at  a deformation 
rate  of  107  s . The  weak  dependence  of  the  viscosity  on  the  nature  of 
the  substance  was  also  reported  in  References  1-3,  which  showed  that,  under 
dynamic  loading  conditions  characterized  by  deformation  rates  of  6 x 10^ 
to  6 x 103  s'1,  aluminum  and  various  grades  of  steels  behave  as  viscous 
liquids  with  y = 3-4  x 104  Pa-s. 

4ir^  a D 

The  deformation  rate  given  by  Mineev21  was  of  the  order  e = y- — , 

8 X 

and  the  viscosity  of  aluminum  at  310  x 10  Pa  remained  constant,  within 
the  limits  of  the  experimental  scatter,  for  deformations  rates  between 

■ In  Mineev 's21  formula  for  the  dynamic  viscosity,  the  initial  amplitude 
a0  does  not  appear . The  a0  was  inserted  by  the  author  to  enable  the 
viscosity  to  have  the  correct  units.  The  term,  k = 2 tt/\,  is  constant 
according  to  Reference  21,  but  it  is  unclear  as  to  which  value  of  X (A.  or 
XoJ  is  to  be  used.  However  x = S(t)/X  and  the  term,  A x/k,  results  In  1 
the  cancellation  of  A,  However  k must  be  nondimensional,  and  it  is 

2m 

assumed  that  Mineev  meant  k to  be  — • 
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4 x 108  to  8 x 10^  s"  21 *23.  This  would  suggest,  from  the  available 
experimental  data,  that  the  viscosity  of  aluminum  is  independent  of  the 
rate  of  deformation  and  hence  the  viscosity  is  Newtonian.  Mineev21  also 
showed  that  aluminum  exists  in  the  plastic  state  behind  the  disturbed 
shock  front . 


Thus,  the  viscosity  of  aluminum  reported  in  Reference  21  is  an  order - 
of -magnitude  smaller  than  the  value  reported  in  References  1 - 3 at  310 
x 10®  Pa  due  solely  to  the  difference  in  temperature  between  the  two 
experiments.  When  the  shock  wave  pressure  is  increased  from  310  x 10® 

Pa  to  1050  x 108  Pa,  the  viscosity  of  aluminum  increases  slightly  but 
remains  less  than  104  Pa-s.  More  exact  values  of  the  viscosity  were  not 
obtained  since  the  Value  of  Ax  could  not  be  resolved  smaller  than  0.02  - 
0.03.  Also,  there  was  evidence  that  large  pressure  gradients  behind 
the  disturbed  shock  front  were  present  due  to  large  initial  disturbances. 

Mineev20  used  the  same  experimental  arrangement  to  obtain  the  dynamic 
viscosity  of  water  and  mercury  under  shock  loading  conditions  at 
pressures  between  40  x 108  and  440  x 108  Pa.  The  viscosity  was  approxi- 
mately the  same  (~  103  Pa-s)  as  the  value  for  shock  compressed 
solids1-3*21*23.  For  both  water  and  mercury,  the  coefficient  of  viscosity 
increases  with  increasing  pressure  (or  increasing  density)  and  decreases 
with  increasing  temperature.  Under  shock  loading  conditions,  the  density 
and  temperature  increase  simultaneously  making  it  difficult  to  accurately 
determine  the  effect  of  only  pressure  or  temperature  on  the  dynamic 
viscosity.  The  viscosity  values  of  water  and  mercury  under  the  test 
conditions  of  Mineev20  are  listed  in  Table  II.  From  the  table,  it 
follows  that  the  deformation  was  of  the  order  20  to  140  percent  and  the 
rate  of  deformation  was  1-6  x 108  s'*.  Under  these  conditions,  in  a 
range  of  pressures  from  80  - 105  x 108  Pa  and  temperatures  of  700  - 1200°K, 
the  coefficient  of  dynamic  viscosity  of  water  is  practically  constant 
(p  ~ 103  Pa-s) . The  viscosity  values  of  water  and  mercury  are  about  the 
same  as  the  viscosity  of  shock-compressed  aluminum  and  lead  at  pressures 
of  300  - 400  x 108  Pa.  It  is  significant  that  the  viscosity  values  for 
aluminum  and  lead  were  obtained  at  approximately  the  same  values  of 
deformation  and  deformation  rate  as  those  of  shock-compressed  water  and 
mercury20*21*23.  According  to  Reference  20, this  explains  the  approximate 
equality  of  the  coefficients  of  dynamic  viscosity  of  shock-compressed 
aluminum,  lead,  sodium  chloride,  water,  mercury,  and  other  materials, 
cited  in  References  21  and  23,  which  are  very  different  under  normal 
conditions.  The  number  of  defects  produced  during  the  deformation  of 
solids  depends  weakly  on  the  type  of  material  undergoing  deformation 
and  is  basically  determined  by  the  deformation  and  deformation  rate. 

Since  the  experiments  of  References  20,  21,  and  23  had  nearly  equal 
values  of  the  deformation  and  deformation  rate,  equal  numbers  of  defects 
can  be  expected  to  appear  and,  as  a result,  approximately  equal  values 
of  the  dynamic  viscosity  coefficients  can  be  expected. 
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TABLE  II.  VISCOSITY  OF  WATER  AND  MERCURY, MINEEV 
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At  a pressure  of  40  x 108  Pa,  the  coefficient  of  dynamic  viscosity  of 
water  is  at  least  an  order-of-magnitude  smaller  than  the  viscosity  value 
at  80  x 108  Pa.  Under  certain  thermodynamic  conditions  behind  the  shock- 
wave  front,  the  development  of  perturbations  does  not  depend  on  X,  and 
this  indicates  a sharp  decrease  in  the  viscosity  of  the  material.  In 
other  words,  i.1  shock  wave  front  breaks  down  due  to  the  increasing 
viscosity  of  the  water  in  the  shock-wave  front.  Also,  at  40  x 10°  Pa  of 
pressure,  the  water  may  exist  in  the  liquid  and  ice  phase,  whereas  the 
water  is  completely  liquid  above  a pressure  of  40  x 108  Pa  20 . 

It  is  also  worthy  of  noting  that  the  viscosity  coefficient  depends 
on  the  shock  relaxation  time,  x ~ or  n = pr  (C02  - C«,  )>  where  C0 
and  C„  are  the  equilibrium  and  non^equi librium  speeds  of  sound.  For 
C 2 » C„2,  p = 1 Mg/m3,  C0  ~ 106  mm/s,  and  t ~ 10~5s,  y ~ 104  Pa-s, 
which  is  in  approximate  agreement  with  the  experimental  results20. 


Mineev20  also  points  out  that  the  shock-wave  thickness  varies 
strongly  with  pressure,  and  the  viscosity  of  water  can  vary  from  10'3 
to  0.3  Pa-s  for  pressures  between  0 to  20  x 10^  Pa. 

Harlow  and  Pracht27  of  LASL  report  a viscosity  coefficient  of  about 
104  Pa-s  for  shock  loaded  aluminum  in  the  later  stages  of  the  experiment, 
after  some  cooling  had  taken  place.  The  viscosity  coefficient  could  drop 
two  orders-of-magnitude  corresponding  to  maximum  heating  of  the  metal. 
Also,  in  Reference  27  the  dynamic  viscosity  of  iron  was  deduced  to  be 
around  2 x 103  Pa-s.  This  value  is  of  the  right  order-of-magnitude  to 
explain  the  decrease  in  jet  velocity  near  the  critical  angle  for  a 
collapsing  shaped-charge  liner27. 

Burkhardt26  deduced  the  dynamic  viscosity  of  liquid  metals  near  the 
melt  temperature  and  found  y = 10-3  * * Pa-s  for  aluminum,  copper,  and 
steel  from  an  investigation  of  the  boundary  surfaces  of  explosively  welded 
metals  at  various  Reynolds  numbers. 

Specifically,  for  liquid  metals  near  the  melting  temperature,  y was 
2.9  x 10'3  Pa-s  for  aluminum,  3.3  x IQ'3  Pa-s  for  copper,  and  2.8  x 10’3 
Pa-s  for  steel, 

Godunov18  investigated  the  viscosity  of  metals  under  impact  in  the 
explosive-welding  regime.  By  means  of  the  method  of  fixed  lines,  it  was 
shown  that  the  viscosity  is  inversely  proportional  to  the  particle 
displacement  in  the  direction  of  the  contact -point  velocity.  The  particle 
displacements  were  measured  by  optical  metallography.  Values  of  the 
viscosity  coefficients  were  estimated  for  aluminum,  copper,  and  steel. 

The  experimental  setup  from  Reference  18  is  shown  schematically  in 
Figure  2a.  Two  plates, 1 and  3, of  the  same  metal  were  prepared  as 
follows,  In  the  upper,  usually  thinner  plate  3,  they  drilled  a hole 
0.3-0. 5 mm  in  diameter  into  which  a wire  4 of  the  same  material  was 
tightly  pressed.  Into  a rectangular  slot  extending  across  the  entire 
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width  of  plats  1,  they  inserted  a plate  2 of  exactly  the  same  thickness 
as  the  slot.  To  improve  the  definition  of  the  boundary  and  eliminate 
air  gaps,  steel  specimens  1 and  2 were  soldered  with  copper  and  copper 
specimens  with  silver  solder  in  a vacuum.  The  space  between  aluminum 
specimens  1 and  2 was  filled  with  epoxy  resin.  The  thickness  of  the 
filler  did  not  exceed  O.OSmm, 


The  plates  i and  3,  thus  prepared,  were  mounted  on  a wooden  base  7. 
Plate  3 was  accelerated  by  the  explosion  products  resulting  from  the 
detonation  of  the  charge  5 (detonator  6)  and  driven  against  plate  1 
(see  Figure  2A).  The  impact  (see  Figure  2B,  detonation  front  6)  usually 
resulted  in  the  explosion  welding  of  plates  1 and  3. 


The  monolithic  specimens  obtained  after  welding  were  cut  in  the 
direction  of  motion  of  the  contact  point  velocity  and  thin  sections  were 
prepared.  On  macrophotographs  of  these  sections,  the  horizontal  displace- 
ment, z,  of  the  boundary  with  respect  to  its  original  position,  was  mea- 
sured as  a function  of  the  distance  to  the  interface  between  the  two 
welded  plates,  y18. 

• The  conservation  of  momentum  for  the  collision  of  two  plates,  the 
experimental  dependence  of  z versus  y,  and  the  steady-state,  incompress- 
ible, viscous,  Navier-Stokes  equation  enables  one  the  estimate  the  metal 
viscosity  coefficients18.  These  viscosity  coefficients  increase  from 
aluminum  through  copper  to  steel.  The  maximum  value  of  the  viscosity 
coefficient  for  steel  (steel  3)  is 

y = (3.9-4. 8)  x 104  Pa-s, 

for  copper  (M3)  the  viscosity  is  about  one-half  as  large: 

y = (2-2.7)  x 104  Pa-s, 

whereas  for  aluminum  (D16)  the  viscosity  is  an  order-of -magnitude  less 
than  for  steel: 


y = (0.31-0,86)  x 104  Pa-s. 

These  viscosity  values  for  aluminum  coincide  with  the  results  of 
References  21  and  23,  and  the  viscosity  values  for  steel  agree  with  the 
results  of  References  2 and  3.  The  viscosity  coefficients  obtained  dis- 
agree with  the  results  of  References  1-3  for  aluminum,  which  predict  a 
value  of  the  order  of  10^  Pa-s,  approximately  the  same  as  for  steel,  and 
there  also  exists  a discrepancy  with  the  data  of  References  21  and  23  for 
steel.  In  References  21  and  23  the  viscosity  of  aluminum,  copper,  and 
steel  are  the  same  (103  Pa-s)  but  an  order- of-magnitude  less  than  the 
values  given  in  References  1-3.  According  to  Godunov18,  the  data  of 
References  21  and  23,  which  found  the  viscosity  values  of  very  different 
materials  to  be  nearly  the  same,  are  suspect.  In  the  opinion  of 
Reference  18,  the  constancy  of  the  flow  behind  the  front  of  the 
sinusoidal  shock  wave  in  Reference  23  was  not  always  realized. 


25 


isiassisi 


Next,  Korsunskii22  used  several  theoretical  methods  to  calculate  the 
dynamic  viscosity  of  liquids  displaying  metallic  properties.  This  data, 
presenting  both  theoretical  and  experimental  molten  metal  viscosity 
values,  are  given  in  Table  IIIA.  Table  I I IB  gives  viscosity  data  for 
molten  sodium  as  a function  of  temeprature. 

Mali1^  quotes  kinematic  viscosity  (y/P)  values  taken  from  References 
18  and  37.  The  kinematic  viscosity  was  given  as  2,5m^/s  for  copper, 
2.5nvVs  for  duraluminum  (D16) , 1.0m^/s  for  aluminum,  0.5m2/s  for  lead, 
and  5.5m2/s  for  steel.  By  assuming  a constant  density,  the  dynamic 
viscosities  become  2.2  x 104  Pa-s  for  copper,  0.7  x 104  Pa-s  for  dur- 
aluminum, 0.3  x 104  Pa-s  for  aluminum,  0.6  x IQ4  Pa-s  for  lead,  and 
4.3  x 104  Pa-s  for  steel. 


Table  IV,  taken  from  Godunov21*,  lists  values  of  \>,  the  kinematic 
viscosity,  and  e for  aluminum  and  steel  obtained  by  analyzing  experimental 
data  from  the  generation  of  planar  flows  formed  by  the  oblique  impact  of 
explosively  driven  plates.  The  dynamic  viscosity  is  also  given  in  Table 
IV  and  was  obtained  by  assuming  incompressible  flow.  The  viscosity 
coefficients  were  determined  based  on  the  assumption  that  decreases  in 
the  flow  rate  in  comparison  to  the  rate  predicted  by  jet  collapse  theory 
(for  small  va'.ues  of  the  collision  angle  and  stagnation  point  velocity) 
occur  because  of  the  action  of  viscous  forces.  Godunov21*  also  gives 
formulae  for  the  variation  of  strain  rate  with  temperature  and  for  the 
Maxwellian  viscosity  as  a function  of  temperatui'e. 

Godunov4  summarized  the  viscosity  measurements  discussed  in  the 
references  cited  above.  From  Godunov4,  the  deformation  of  a material 
depends  essentially  on  the  load  acting  on  it  as  well  as  the  rate  of  its 
application.  At  present,  the  range  of  low  strain  rates,  achievable  from 
existing  tension  test  machines,  have  been  investigated  in  detail  where 
the  strain  rate  can  vary  from  10"4  to  103  s”l*  Shock  tests  with  ex- 
plosives permit  strain  rates  of  103  to  10§  s~*  to  be  obtained. 


The  viscoplastic  solid  model  was  used  by  References  1 and  3 in  the 
analysis  of  experiments  designed  to  determine  the  coefficient  of  viscosity 
of  aluminum  and  steel  for  e = IQ3  s . The  method  of  determining  the 
viscosity  of  metals  during  explosive  welding  is  described  in  Reference 
18.  The  strain  rates  used  in  these  tests  were  between  5 x 103  - 3 x IQ4 
s-1. 


In  Reference  5,  a method  was  presented  to  determine  the  viscosity 
of  metals  during  the  symmetric  oblique  collision  of  plates  in  the  mode 
where  reverse  jets  exist.  In  this  case,  the  jet  velocities  obtained 
experimentally  were  compared  with  the  velocities  computed  by  the  viscous 
deceleration  model.  In  this  case,  estimates  of  the  coefficients  of 
viscosity  of  diverse  metals  were  carried  out  in  the  103  - 10'’  s”1  range 
of  strain  rates.  The  results  obtained  under  the  conditions  of  oblique 
collisions  of  metal  plates  are  in  good  agreement  with  the  results  of 
References  21  and  23,  where  the  development  of  small  perturbations  on 
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TABLE  I I I A.  CALCULATED  AND  EXPERIMENTAL  VISCOSITIES 
OF  SEVERAL  MOLTEN  METALS  NEAR  THE  MELTING  POINT. 
(Reference  22) 


Metal 


y Theoretical  x 10°  Pa-s  y Experimental  x 10  Pa-s: 


Li 

.28  - .44 

.59 

K 

.83  - 3.6 

.52 

Mg 

2.0 

1.3 

Zn 

1.1  - 2.9 

2.8 

Cd 

3.6 

2.4 

Hg 

1.64 

1.54 

A1 

1.8  - 2.5 

1.2 

Ga 

1.9 

1.9 

In 

4.3 

1.65 

Si 

2.9 

2.0 

Sn 

1.7 

1.8 

Pb 

4.8  - 6.1 

2.6 

Sb 

3.8 

1.5 

Bi 

2.4 

1.7 

TABLE  IIIB.  VISCOSITY  OF  MOLTEN  SODIUM  AS  A 
FUNCTION  OF  TEMPERATURE. (Reference  22) 

3 3 

Temperature  (°K)  y Theoretical  x 10  Pa-s  y Experimental  x 10  Pa-s 


373 

1.4 

- 2.7 

0.69 

443 

1.1 

- 2.0 

0.5 

473 

0.9 

- 1.66 

0.44 

513 

0.8 

- 1.5 

0.39 
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the  shock  front  was  studied  for  the  same  range  of  strain  rates. 

Figure  3,  taken  from  Godunov4,  shows  the  variation  of  experimental 
viscosity  coefficients  of  steel  and  aluminum  with  strain  rate.  For 
strain  rates  greater  than  104  s~^ , the  viscosity  coefficient  of  aluminum 
and  steel  depends  weakly  on  strain  rate,  i.e.,  the  viscosity  deceases 
by  about  a factor  of  ten  as  the  strain  rate  increases  from  10 • s 1 to 
lO?  s'1.  Thus,  the  metals  can  be  considered  approximately  as  Newtonian 
fluids.  The  viscosity  increases  rapidly  as  the  strain  rate  decreases 
below  104  s . The  viscosity  data  in  Figure  3 was  taken  from  References 
1,  3,21  and  23  at  strain  rates  above  10*  s.  The  viscosity  values  at 
strain  rates  below  10^  s'1  were  taken  from  Reference  39  for  tests  in- 
volving Tolling  of  the.  metal  with  preheating  to  1273°K  and  stamping  of 
the  metal  at  1273°K  as  reported  in  Reference  4. 

Other  USSR  viscosity  data  was  reported  by  Ilakhsiyan1^,  who  quoted  a 
dynamic  viscosity  coefficient  of  armor  steel  as  3.9  x 104  Pa-s.  Also, 
Ivanov25  studied  the  plastic  deformation  of  soft  steel  tubes  by  detonat- 
ing spherical  charges  of  explosive  placed  in  the  central  cross  section 
of  the  tubes.  The  strain  rates  obtained  in. this  test  ranged  from  10->  to 
lO?  s'1.  The  viscosity  values  obtained  by  Ivanov25  were  claimed  to  be 
in  agreement  with  the  viscosity  values  obtained  by  other  methods  such 
as  References  1,  2 and  57.  The  soft  steel  viscosity  values  were  given 
as  2-5  x 104  Pa-s  by  Reference  25+. 

Shlykov-8  found  that  the  kinematic  viscosity  of  copper  alloys  varied 
strongly  with  temperature  and  composition.  These  copper  alloys  were 
investigated  experimentally  in  the  liquid  stage.  The  alloys  employed  were 
Cu-Zn-Si,  Cu.-Zn,  and  Cu-Si  with  up  to  50%  Zn  and  up  to  12%  Si.  The 
concentration  dependence  of  Si  (in  Cu-Si  alloys)  on  the  kinematic 
viscosity  is  quite  complex.  Small  concentrations  of  Si  increase  the 
kinematic  viscosity.  Larger  amounts  of  Si  cause  the  viscosity  to  de- 
crease, and  still  larger  amounts  of  Si  cause  the  kinematic  viscosity  to 
increase  again38. 

The  US  investigators:  LASL,  SANDIA  Laboratory,  and  Lawrence  Livermore 
Laboratory,  have  measured  dynamic  viscosities  under  explosive  loading 
conditions  to  be  of  the  order  102  Pa-s,  typically  lower  than  the  mea- 
sured USSR  viscosity  values.  Swanson28  studied  the  effects  of  strain- 
hardening  and  strain-rate  treatments  on  one-dimensional  calculations  of 
aluminum  response  during  plate  impact  experiments.  For  strain-rate 
effects  on  aluminum,  a dynamic  viscosity  of  200  Pa-s  was  shown  to  best 
represent  the  free-surface  velocity  versus  time  behavior  and  the 
observed  attenuation  of  thin  pulses  by  rarefaction  waves  during  the 


f The  units  given  by  Ivanov2'0  are  kg  f/m 2,  which  are  not  dynamic  viscosity 
units.  The  author  has  assumed  that  Ivanov  means  kg  f-s/my  Viscosities 
of  order  10 4 Pa-s  are  then  in  agreement  with  the  viscosities  given  by 
other  USSR  scientists  as  claimed  by  Ivanov 25 . 
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plate  impact  experiments.  Strain-hardening  effects  were  shown  to  be 
secondary,  but  complementary,  to  the  effects  of  viscosity.  Pressure  and 
temperature  dependent  relationships  for  the  yield  strength  and  shear 
modulus  showed  only  minor  effects  and  were  concluded  to  be  of  secondary 
importance.  Swanson29  noted  that  the  dynamic  viscosity  value  is  asso- 
ciated with  the  movement  of  dislocations  that  occur  only  after  the  yield 
stress  has  been  reached.  Dislocation  dynamics,  for  shock  impact  experi- 
ments, have  been  studied  by  many  USA  researchers,  notably,  Asay40*41, 
Taylor42,  Gilman43*44,  Smith45,  and  Lipkin46.  In  the  USSR  crack  propaga- 
tion and  dislocation  theory  was  studied  by  Finkel47*48  and  Golovin49. 

; Giiraan30  theoretically  calculated  the  approximate  dynamic  viscosity 
for  strong  shock  waves  in  metals  and  concluded  that  the  dynamic  viscosity 
should ;be  of  the  order  of  0.1  Pa-s,  whereas  the  viscosity  of  liquid  metals 
is  of  the  order  of  5 x 10“^  Pa-s. 

Finally,  Barnes28  in  Figure  4,  shows  the  variation  of  dynamic  viscos- 
ity with  pressure.  Barnes  comments  that  viscosity  would  appear  to  be 
important  in  mitigating  shock  perturbation  growth,  but  the  value  of  the 
-viscosity  coefficient  is  highly  uncertain.  In  Figure  4,  the  viscosity 
data  at  high  pressures  is  from  the  USSR  experiments21 , and  the  data  at 
200  Pa-s  is  from  Swanson29.  The  point  below  Swanson' s results  from 
fitting  Barnes' experimental  data  on  Taylor's  growth  in  aluminum50  with 
a viscous  flow  model.  G.  N.  White  actually  used  the  viscous  flow  model 
to  fit  Barnes' data  (see  Reference  28).  Also,  recall  that  Gilman30 
calculated  an  upper  limit  of  the  viscosity  in  stressed  metals  to  be 
0.1  Pa-s. 

Asay51  has  deduced  stressed  metal  dynamic  viscosities  by  measuring 
shock  rise  times  via  velocity  interferometry.  The  viscosity  values 
given  below  represent  an  upper  limit  based  on  shock  transition  times." 

The  time  resolution  is  limited  to  1 to  3 nanoseconds.  The  maximum  stress 
is  assumed  proportional  to  the  strain  rate,  o ~ u e,  and  the  maximum 
stress  is  determined  between  the  Rayleigh  line  and  the  shock  Hugoniot 
pressure.  The  strain  rate  is  calculated  as  the  strain  induced  by  the 
shock  divided  by  the  shock  rise  time,  which  is  resolution  limited.  The 
calculated  viscosity  values  are  listed  below. 


Material 

Cu 

A1 

4340  Steel 


Pressure,  Pa  x 10 

- 930 
~ 450 

- 930 


-8 


Dynamic  Viscosity,  Pa-s 

- 3 x 102 

- 1 x 102 . 

~ 6 x TO2 


Also,  Harrison52  is  currently  computer  modeling  the  USSR  high-pres- 
sure experiments  in  order  to  resolve  the  strong  variation  in  the  quoted 
viscosity  values  especially  regarding  the  value  of  104  Pa-s  quoted  by 
the  USSR  versus  the  IQ2  Pa-s  viscosity  values  quoted  by  the  USA. 
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VISCOSITY 


Additional  data  in  the  open  literature  considers  factors  related  to 
viscosity,  although  viscosity  coefficients  are  not  directly  employed. 
These  factors  are  thermal  effects  during  impact  from  Belyakov53* 5**,  the 
effect  of  temperature  on  the  mechanical  properties  of  metals  from 
Sazonov55  and  Sklyuev56,  explosive  effects  on  plate  acceleration  from 
Godunov57  and  Trofimov55,  and  the  mechanical  properties  of  metals  under 
dynamic  conditions  from  Yakhnin5^,  Tass66,  Men'shikov51, and  Wilkins62. 
References  61  and  62  treated  elastic-plastic,  stress -strain  models  with 
strain  hardening  considerations  in  order  to  determine  the  dynamic  yield 
strength  of  rods  of  various  materials  impacting  rigid  targets.  Wilkins 
lists  the  calculated  dynamic  yield  strength  of  several  materials  in- 
cluding aluminum,  steel,  copper,  uranium,  tantalum,  magnesium  alloys, 
and  beryllium  alloys. 


It  is  interesting  to  note  that  Moss63  of  BRL/ARRADCOM  experimentally 
obtained  viscosity  data  from  void  growth  measurements.  Dr.  Moss 
obtained  dynamic  viscosity  values  of  357  Pa-s  for  1/2-inch  thick  (12,7 
mm)  RHA  plate  and  20  Pa-s  for  copper.  The  viscosity  value  for  the  RHA 
plate  includes  some  effects  of  void  coalescence  but  is  in  general 
agreement  with  the  value  given  by  Asay51.  The  viscosity  value  for 
copper  is  lower  by  about  an  order-of-magnitude  over  Asay's  value. 
However,  the  exact  test  conditions  (thermodynamic  properties)  of  the 
test  specimens  are  not  known. 


IV.  SUMMARY  AND  CONCLUSIONS 


The  USA,  for  the  most  part,  assumes  elastic-perfectly  plastic  or 
strain  hardening  models  for  the  constitutive  relationships  used  in  the 
hydrocode  calculations.  The  USSR,  on  the  other  hand,  contends  that 
viscous  effects  are  at  least  as  important  as  compressibility  effects 
in  the  jet  collapse  process.  Thus,  the  USSR  constitutive  relationships 
utilize  a visco-plastic,  stress-strain  rate  model  where  the  stress  is 
assumed  proportional  to  the  strain  rate  and  with  the  factor  of 
proportionality  being  the  dynamic  viscosity  coefficient.  The  viscosity 
coefficient  is  assumed  to  be  constant,  but  this  is  true  only  fcr  a given 
pressure,  temperature,  and  strain  rate. 


The  USSR  measured  the  viscosity  coefficients  for  the  propagation 
of  shock  waves  through  several  materials.  The  viscosity  of  all  shock 
loaded  materials  is  of  the  same  order-of-magnitude.  For  metals,  the 
dynamic  viscosity  coefficient  was  in  the  range  of  10^  to  10$  Pa-s. 

The  US  investigators  typically  measure  dynamic  viscosities  of  the 
order  of  10^  Pa-:;.  The  estimate  of  the  viscosity  coefficient  from 
measured  values  undoubtedly  depends  on  the  experimental  procedure 
employed  as  well  as  the  experimental  accuracy.  Theoretical 
viscosity  coefficients  are  about  0.1  Pa-s,  and  liquid  metal  viscosity 
coefficients  are  of  the  order  10-3  Pa-s.  The  differences  in  the  viscos- 
ity values  between  the  USSR  and  USA  might  be  attributed  to  pressure 
differences,  Barnes28  (see  Figure  4).  Strain  rate  also  influences  the 
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viscosity  as  shown  by  Godunov1*  (see  Figure  3) . Temperature  also 
inf luences the  viscosity  strongly,  especially  near  the  melt  point. 

In  short,  as  stated  by  Barnes28,  "Viscosity  would  appear  to  be  a 
very  important  factor  in  mitigating  perturbation  growth,  but  the  value 
of  viscosity  coefficient  is  highly  uncertain." 

It  seems  to  this  author  that  further  research  in  this  area  is 
required  to  determine  the  dependence  of  the  viscosity  coefficient  on 
pressure,  temperature,  and  strain  rate  as  well  as  to  increase  the 
accuracy  of  the  viscosity  measurement.  Viscous  effects  should 
definitely  be  an  important  mechanism  in  dissipating  shock -wave  fronts 
during  the  shaped-charge  jet  formation  and  the  jet  penetration  process. 

Viscous  effects  can  be  incorporated  into  the  jet  calculations  (for 
example,  Walters64),  but  the  coefficient  of  dynamic  viscosity  must  be 
known.  The  problem  is  further  complicated  by  the  dependence  of  the 
viscosity  on  both  the  pressure  and  temperature,  since  good  equation-of- 
state  data  and  accurate  temperature  calculations  are  currently  not  avail- 
able for  problems  involving  the  explosive  loading  of  metals. 


V.  RECOMMENDATIONS 

Visco-plastic  or  rate- dependent  models  should  be  investigated  in 
regard  to  the  possibility  of  improving  shaped-charge  jet  collapse, 
formation, and  penetration  calculations. 

The  one-dimensional,  incompressible,  steady-state  model  for  the  jet 
velocity  given  by  Godunov4,  which  is  analogous  to  the  Pugh-Eichelberger- 
Rostoker  model  except  that  viscous  effects  have  been  included,  should  be 
analyzed.  In  particular,  the  jet-no-jet  criterion  proposed  by  Godunov, 
which  is  based  on  a critical  Reynolds  number,  should  be  checked  against 
the  USA  data.  A critical  Reynolds  number  criterion  on  the  jet-no- jet 
conditions  is  more  meaningful  than  a critical  Mach  number  criterion  for 
incompressible  flow  models.  The  necessary  equations  are  given  in  this 
report . 
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